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We present an isospin analysis of the decay modes B — ► DD, D*D, DD* and D*D* , 
^h ; 
Qh which allows determination of the final-state interaction phases. As these transitions 

Mh. have branching ratios of the order 10~ 4 or larger, they could be useful in detecting 



i-G the CP-violating phase 4>\ = arg(— V^VcdVtbV^) in the first-round experiments of a 

i?-meson factory. The problem of penguin pollution may still be present. Once the 
Kobayashi-Maskawa matrix elements are known, it is possible to obtain the magnitudes 
and relative phases of hadronic matrix elements for B — ► D^*>D^*> . This will in turn 
lead to some information about the penguin pollution. 
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1 Introduction 

Within the standard electroweak model, three angles of the Kobayashi-Maskawa (KM) uni- 
tarity triangle 

v: b v ud + v; b v cd + v t iv td = o, (i) 

denoted by fa (i = 1, 2, 3), are determined from CP asymmetries in neutral 5-meson decays 
to hadronic CP eigenstates |l|, [|, ||. It is expected that fa = a.rg(—V* b V cd VtbV td ) can be 
unambiguously extracted from the decay rate difference between B® — > ipKs and B d — ► ipKs- 
While this decay mode should be sufficient for determining 0i, the initial luminosity of a 
B factory may require us to search for some additional decay channels which could help us 
establish the presence of CP violation as quickly as possible Q. For this purpose, we shall 
investigate 

B — ► DD , D*D , DD* , D*D* (2) 

in some detail. 

In practical experiments the decay mode B d — > D + D~ should have fairly large branching 
ratio. Under SU(3) symmetry, one can make the following rough estimation: 

B(B° d -> D + D~) ~ sin 2 C B{B° d -► D+D~) « (3.4 ± 1.9) x 1(T 4 , (3) 

where 9q is the Cabibbo angle, and B(B d —>■ DfD~) = (7 ± 4) x 10 -3 has been measured 
in experiments ||. Also, the penguin effect in B d — » D + D~ is expected to be smaller than 
that in B d — > 7t + tt~ . Thus the CP asymmetry between B® — > D + D~ and B® — * D + D~ may 
be dominated by angle 4>\. In contrast with B d — > D + D~ , the decay modes S^ — > D + D*~ , 
D* + D~ and D* + D*~ undergo the same weak interactions and have the comparable branching 
ratios, although they are not the exact CP-even eigenstates ||. 

Of course the measurement of CP violation in B d — > D + D~ can not only cross-check 
the extraction of <p\ from B d -^ ipK$, but also shed some light on the penguin effects and 
final-state interactions (FSIs) in nonleptonic B decays to double charmed mesons. For this 
reason, it is worth studying both B d — > D + D~ and B d — > D°D° in a model-independent 
approach. The similar treatment is applicable to the processes B d — ► DD*, D*D, etc. 

In this work we shall carry out an isospin analysis of the processes B — > D^D^*\ to relate 
their weak and strong phases to the relevant observables. It is found that the time-dependent 
measurements of B d —* D + D~ and B d —>■ D°D° together with the time-independent mea- 
surements of B£ —> D + D° and B~ — > D~D° allow one to extract a phase parameter <f)[, 
which consists of both <fii and the penguin-induced phase information. Direct CP asymme- 
tries in B d — > D + D~ and D°D° are time-independently detectable on the T(4S) resonance. 



For numerical illustration, we apply the effective weak Hamiltonian and factorization ap- 
proximation to B+ — ► Z)M + Z)(*)° and B~ — ► /jM^/jM , since each of them is only involved 
in a single isospin amplitude. We find that their branching ratios are all above 10~ 4 and the 
relevant time-independent CP asymmetries may reach the 3% level. The time-dependent 
CP asymmetries in B d — ► £)(*)+Z)(*) - and £)W°D(*)° are expected to be of order 1. We 
also emphasize that it is possible to obtain some information on the magnitudes and relative 
phases of hadronic matrix elements in the isospin approach, once the KM matrix elements 
have been known. 

2 Isospin analysis 

The effective weak Hamiltonians, responsible for B~ — > D~D°, B d — »■ D + D~, B d — »■ D°D° 
and their CP-conjugate processes, have the isospin structures 1 1/2, — 1/2) and 1 1/2, +1/2) 
respectively. The decay amplitudes of these transitions can be written in terms of the isospin 
amplitudes: 
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A+- = (D + D-\7i eS \B° d ) 


= 2^i + Ao 


A 00 = (D°D \n cS \B° d } -. 


- \ (Ai - A ) 


A +0 = (D + D \H cS \B+) 


= A i; 



(4) 



and 



(D + D-\H eS \B° d ) 


= \ {Ai + A 


(D D°\n cS \B° d ) -- 


= \(A,- A 


(d-d°\h cS \b-) 


= A,. 



A + - E 

A 00 E 

A-° = (D-D°\H cS \B-) = A, . (5) 

Here A\ {A\) and Aq (Aq) are the isospin amplitudes with 1=1 and 1 = 0, respectively. 
Clearly the isospin relations (4) and (5) can be expressed as two triangles in the complex 
plane (see Fig. 1 for illustration): 



roo _ T-0 



A + ~ + A w = A" . 



A + - + A w = A- u . (6) 

One is able to determine the relative size and phase difference of isospin amplitudes A\ (A\) 
and Aq (Aq) from the above triangular relations. Denoting 

Ao _ ie Aq _ q 

A = ze ' A = ze ' (7) 



A* 





A-° 
(a) (b) 

Figure 1: The isospin triangles of B — > DD in the complex plane. 



then we obtain 
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and 
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2{\A+-\ 2 + \A 00 \' 
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arccos 



\A 



0012 s 
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If z = 1 and 9 = 0, for example, we find that \A 00 \ = 0, i.e., the decay mode B c 
forbidden. 



(9) 



D°D° is 



Note that 9 (9) is in general a mixture of the weak and strong phase shifts, since both 
Aq (Aq) and A\ (A\) may contain the tree-level and penguin contributions. This point can 
be seen more clearly if one writes the isospin amplitudes Ai and Aj (I = 1, 0) with the help 
of the low-energy effective AB = ±1 Hamiltonians. For example, Aj can be given as 

G, 



with 



Aj = ((DDMTUaiAB = +1)\B) = -| E [(V* b V qd )Sj 
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Sj = ci((D5)j|Q?|S> + c 2 ((DD)j\Q q 2 \B) + ^ [c^D^AQ^B) 



(10) 



(11) 



i=3 



where Wilson coefficients q and four-quark operators Qi at the scale \x = 0{rrib) have been 
well defined in Ref. 0. The expression of Aj is straightforwardly obtainable from Eq. (10) 
through the replacement (V* b V q d) — > (V q bV* d ). The tree- and penguin-type hadronic matrix 
elements in Sf are expected to consist of different strong phases, and these phases should be 
different from those in Sf. This implies that the overall phases of A\ \A\) and A (A ) are 
non-linear combinations of the same weak phases and the different strong phases, therefore 
9 (9) is neither purely weak nor purely strong. 



Finally it is worth mentioning that the same isospin relations hold for the decay modes 
B — > DD* and B — > D*D. Of course, the isospin parameters z (z) and (0) in B — ► 
1)5, DD* and £)*£) may be different from one another due to their different FSIs. As for 
B — ► D*D*, the same isospin relations hold separately for the decay amplitudes with helicity 
A = -1, 0, or +1. 

3 Time-independent measurements 

The quantities |>1 +0 | and \A~°\ are obtainable from the time- independent measurements 
of decay rates of B+ -> D + D° and B~ -> D~D°. A determination of \A + ~\ (\A 00 \) and 
|A +_ | (|^4 00 |) is possible through the time-integrated measurements of B d vs B d — > D + D~ 
(D°D°) on the T(4S) resonance, where the produced two B d mesons are in a coherent state 
(with odd charge-conjugation parity) until one of them decays. In practice, one can use the 
semileptonic transition of one B d meson to tag the flavor of the other meson decaying to 
D + D~ or D°D°. The probability for observing such a joint decay event reads H. ffl: 



\A+-\ 2 + \A+-\ 2 __ 1 \A+-\ 2 -\A + -^ 
1+^ 



K^X^D+D-) oc |^| 2 [J LZJ L T __ . J L_J !_ , (12) 



or 

1±YT. nOnO\ _ I A 12 



MOO 1 2 _|_ MOO 1 2 | M 00 | 2 — l/^ 00 ! 2 ' 

~2 T l + x 2 d 2 



TZ^X^D^D") oc |A| 2 J — ~ir — ~ T ^^■ J — Sr — L , (13) 



where \A t \ = \(l + X-\H e g\B%)\ = \ (l~ X + \H cS \ B%) | under CPT symmetry, andx d = Am/T « 
0.73 is a measure of S^ — -B° mixing ||. By now the semileptonic B d transitions such 
as B d — > D^*>~V r v l and B^ — > D^*' + l~v l have been well reconstructed ||, i.e., \A{\ has 
been detected independent of the above joint decay modes. Once 1Z{l ± X^\D + D~) and 
7?.(/ ± X =F ; D°D°) are measured, we shall be able to determine the quantities \A^ | (|^4 00 |) 
and \A+-\ (\A 00 \). 

The time-independent measurements mentioned above allow one to construct the isospin 
triangles in Fig. 1. Consequently the isospin parameters z (z) and 9 (6) are extractable in 
the absence of any time-dependent measurement. If the branching ratios of B d — > D°D° 
and B d — > D°D° are too small to be observable, then large cancellation between the isospin 
amplitudes A\ [A\] and A (A ) must take place. In the case that B d — > D + D~ and 
B£ — » D + D° have been measured earlier than B d — > D°D°, a lower bound on the rate of 
the latter decay mode is model-independently achievable from the isospin relations obtained 
above. Since cos6 < 1, we get from Eq. (8) that 
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B(B° d ^D°D°) > 



B(B° d -> D+D-) _ 



B(B+ -> D + D») , (14) 



where tiny isospin-violating effects induced by the mass difference m D0 — %- and the 
life time difference t b — r B have been neglected. This bound should be useful to set a 
limit for the results of B(B d — »■ D°D°) obtained from specific models of hadronic matrix 
elements. Following the same way, one can find the lower bounds for the branching ratios of 
B° d -»■ D*°D°, D°D*° and D*°D*°. 

The nonvanishing asymmetry between the decay rates of B* — ► D + D° and B~ — ► D~D° 
signifies the existence of direct CP violation. By use of the isospin amplitudes in Eqs. (10) 
and (11), we obtain the following CP asymmetry: 

A±0 = K(B+ -+ D+EP) + n{B- -+ D~D°) ~ ^^ N n ' (15) 

where 03 = arg(— V* b V u dV c bV* d ) is an angle of the KM unitarity triangle, and Nu can be 
read from 

Ntj = KRe(S?SJ*) + kT 1 Re^Sf ) - cos0 3 Re(S?S]* + S?S?) (16) 

with k = \V u bV u d\/\V c bV c d\. For the processes Bd — »■ D + D~ and D°D°, pure signals of direct 
CP asymmetries may manifest themselves on the T(4S f ) resonance: 

TZ(l-X+;D+D-) - K(l + X-;D+D-) 



A, 



and 



A 



00 



2 sin 03 


;D+D~) + TZ(l+X-;D+D-) ' 
Im(5^5'f + S2f5g* + Sj'Sg* + SftSf ) 




JVn + A^oo + A^io + N 01 
■ ) D°D°) - TZ(l + X-;D°D°) 


n(i-x+ 

2 sin 03 


;D°D ) + n(l+X-;D°D°) 

Tm/'Q U C C * ^ QU QC* QU QC* QU QC*\ 



;i7) 



l + xg iVn + N 00 - N w - N i 

If the decay modes B d — »■ D°D° and S° — >■ D°D° are forbidden due to the absence of final- 
state rescattering (i.e., ~ and 2 w 1, or S^ ps So), then measuring the CP asymmetry 
vAoo is practically impossible. In this case, we arrive at an interesting relation between the 
asymmetries A±o and A+-. 

A± « (l+x 2 d )A + - « 1.5.A+-. (19) 

The validity of this relation is testable in the forthcoming experiments at a 5-meson factory. 

It is worthwhile at this point to give a brief comparison between the isospin language 
and the intuitive quark-diagram description for B — ► D^D^*\ Both the isospin amplitudes 
A\ and Ao are dominated by the spectator (external ^-emission) quark graph with the KM 
factor V* b V c d, but they also receive some small contributions from the loop-induced penguin 



and annihilation-type tree quark diagrams. Hence the branching ratios of B% 



£)(*)+£)(*)- 

and B+ — > £)(*)+DM may be of the same order. In the assumption of no final-state 
rescattering or channel mixing, B® — > D^*'°D^*'° takes palce only through the annihilation- 
type quark graphs, which are expected to have significant formfactor suppression in the 
factorization approximation. This argument is compatible with the isospin analysis, since 
the cancellation between A\ and Aq in A 00 implies that the dominant spectator diagram 
does not contribute to B® — > £)(*)°D(*)°. However, one should keep in mind that FSI 
effects are possible to significantly enhance the decay rate of B d — > £)(*)°£)(*) to the level 
comparable with that of B° d — »■ £)(*)+D(*) - or B+ — > £)(*)+ /)(*)° ; making the naive quark- 
diagram language in failure. 



4 Time-dependent measurements 



To probe the CP asymmetry induced by the interplay of direct decay and B® — B® mixing 
in Bd — > DD, the time-dependent measurements are necessary on the T(4S) resonance at 
asymmetric B factories. In such an experimental scenario, the joint decay rates can be given 
as follows [Bl: 



U{l ± X^,D + D--t) oc |^|V r|t| 



\A+-\ 2 + \A + -\ 2 \A+-\ 2 - \A + -\ 2 , 

1 -7T ~ =F J L -= J L cos(x d Tt) 



±\A + -\ 2 Im[± 



qA A 



pA + ~ 



sm(xdTt) 



(20) 



and 



K{l ± X z ? ) D°D°-t) oc |^| 2 e" r|41 

± U 00 | 2 Im 



A 00 1 2 _|_ MOO 1 2 M 00 ! 2 — I A 00 1 2 

— =F cos(xrfrt) 
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'qA^ 
pA 00 



sin(xdTt) 



(21) 



where t is the proper time difference between the semileptonic and nonleptonic decays [], and 

(22) 



q 
v 



expl^-i^oj 



exp(-i2</> J 



stands for the phase information from B® — B d mixing Q. For simplicity, we denote the 
phase difference between Ai and A\ as 



H> 



1 (A, 
- arg — — 

2 S Ui 



1 



arg 



V cb V c * d S?-«S?exp(-i< 



(23) 



V; b V cd Sl-nSfexpi+ifo) 

'•Note that the proper time sum of the semileptonic and nonleptonic decays has been integrated out, since 
it will not be measured at any B-meson factory. 



In terms of the isospin parameters, coefficients of the sin(xdTt) term in Eqs. (20) and (21) 
are given by 



'qA+-\ \A +0 A-^ 



M*— ) = ^W[- Sin (2fi) - zmn(9 + 



and 



+ z sin (6 - 2^;) + zz sin (6-9- 20;)] (24) 



/«4 00 \ l4+°4-°l 

Im ( f A^J = WP [_Sin(20 ' l) + ^ Sin ^ + 2 ^ 



zsin 



in(0-20i) + zzsin (0-0-2$)] , (25) 



where </>; = 0o — <£>. All the quantities on the right-hand side of Eq. (24) or (25), except </>;, 
can be determined through the time-independent measurements of B — > DD on the T(45) 
resonance. Thus measuring the CP-violating observable on the left-hand side of Eq. (24) or 
(25) will allow a model-independent extraction of </>;. 

Let us make two remarks about the results obtained above: 

(1) Within the standard model, 0o ~ <p\ holds to an excellent degree of accuracy. If the 
tree-level quark transition b — ► (cc)d dominates the decay amplitude of B* — > D + D° (i.e., 
\Si\ 3> \Si\ in Eq. (23)), then we get ip ~ a.rg(V c bV* d ) ~ as a pure weak phase. In this case, 
the magnitude of <\>\ is extractable from the time-dependent measurement of Bd — > D + D~ 
or Bd — > D°D° |], Q. Note that a model-dependent estimation in the standard model gives 
ip ~ —3° (see the appendix). It is worth pointing out that the B® — B® mixing phase O 
can be reliably determined from the CP asymmetry in B® vs B^ —>■ ipKs either within or 
beyond the standard model H. Thus a comparison of 0o (extracted from Bd — > ipKs) with 
0; (extracted from Bd — > D + D~ or Bd — ► D°D°) will constrain ip, which may reflect the 
penguin- induced phase information in B — > DD. 

(2) It is interesting to note that we can, in principle, obtain phases and magnitudes of 

the hadronic matrix elements Sj. On the experimental side, \A \ (\A Q \), \Ai\ (|Al|) and 9 

(9) can be determined from the time-independent measurements; and </>; can be extracted 

from the time-dependent measurements. If the KM phases are known, the hadronic matrix 

elements Sq, S™, Sq and Si represent seven unknown parameters as the overall phases of 

them are physically irrelevant. Thus the magnitudes and relative phases of these quantities 

should be determinable from experimental measurements of the relevant branching ratios 

and CP asymmetries. 

§ This point relies on the condition that no significant penguin effect exists in Bd — > i^Ks, and it is 
testable through the time-independent measurement of the decay rate difference between B^ — > ipKs and 
B® — > ipKs on the T(45) resonance. 



A special but interesting case is z — z — 1. It can be obtained if the decay modes 
B — ► DD occur dominantly through the tree- level subprocess b — ► (cc)d or 6 — > (cc)rf. In 
this case, A (A ) and Ax (Ai) have a common KM factor; thus 8 (8) is a pure strong phase 
shift. This will lead, for arbitrary values of 8 and 8, to the results 
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sin - . 
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(26) 



As a straightforward consequence, one gets 

\ A +-\ 2 + |A 00 | 2 
|A+-| 2 + |A 00 | 2 



\A 
IA 



-0|2 



-0|2 



(27) 



i.e., the two isospin triangles in Fig. 1 become right-angled triangles. If 8 = 8 is further 
assumed, we obtain 

'qA + -\ 



Im 



pA H 



T fqA 00 " 
1m — - 



|A+°A-°| , /x 9 ^ 
1 ' sin (200 cos 2 - 
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., sin (20i) sin 2 - . 



(28 



v pA°°; |A°°p 

One can see that these two CP-violating quantities have the quasi-seesaw dependence on 
the isospin phase shift 8. The magnitude of sin^c^) turns out to be 

'gA otr 



sin(2C 



1 



|A+°A-°| 
apparently independent of 8. 



A 



+- 12 



Im 



'gA 



A 00 | 2 Im 



pA 00 



(29) 



5 Summary and conclusion 



We have presented an isospin analysis of the weak decays B 
can be summarized as follows: 



£)(*)£)(*). The main results 



(a) The time-independent measurements of these transitions on the T(4S) resonance 
allow one to extract the isospin quantities and probe the direct CP asymmetries in them. 
It is possible to extract a phase parameter, which consists of the phase information from 
both Bj — Bj mixing and penguin diagrams, from the time-dependent measurements of 
Bd — >■ D + D~ and D°D°. A comparison of this phase with that extracted from Bd — > ipKs 
will be interesting, since their difference signifies the penguin-induced phase information (no 
matter whether new physics is present or not). 



(b) Once the KM matrix elements have been determined, the relevant hadronic matrix 
elements (including their phase information) can be determined, through the isospin analysis, 
from some measurements of the decay rates and CP asymmetries. 

In the appendix, we have made use of the effective weak Hamiltonian and naive fac- 
torization approximation to estimate the branching ratios of B+ — ► £)(*)+£)(*) and B~ — ► 
£)(*)-£)(*)o as we u as their CP asymmetries. It is remarkable that all these decay modes 
can well be detected in the first-round experiments of a S-meson factory. In particular, only 
about 10 8 B^ events are expected to need for the exploration of direct CP-violating signals 
in them (at the 3% level). 

We conclude that a careful experimental study of the decay modes B — ► £)(*)!)(*) at 
the forthcoming B factories, before the measurements of B — ► tttt and other charmless B 
decays become available, will be able to cross-check the extraction of the weak angle 0i from 
Bd —* ipKs, to shed some light on the penguin and FSI effects in B decays to double charmed 
mesons, and to probe direct CP violation in both charged and neutral 5-meson systems. 
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Appendix 

Here we calculate the branching ratios of 5+ — ► £)(*)+/)W and B~ — ► D(*) - D(*)° as well 
as their CP asymmetries numerically, in order to give one a feeling of ballpark numbers to be 
expected within the standard model. It is suitable to apply the effective weak Hamiltonian 
and factorization approximation to these decay modes, because each of them only has a 
single isospin amplitude. In contrast, Bd — ► £)(*)+.D(*) _ or Bd — *• £)W°L)(*) is involved in 
two different isospin amplitudes; thus a direct application of the factorization approximation 
to such transitions may be problematic unless the FSI effects are negligibly small. 

In estimating the branching ratios of 5+ — ► D + D°, D* + D°, D + D*° and D* + D*°, it is 
instructive to neglect small contributions from the hadronic matrix elements (D + D°\Q'i 2 \B+) 
(annihilation) and (D + D°\Q^^io\B^ ) (penguin). These transitions have the weak interaction 
similar to that in £?+ — ► D+D°, D* + D°, D+D*° and D* + D*°, whose decay rates have already 



been measured in experiments ||. Then a comparison between the above two sets of decay 
modes, with the help of the factorization approximation, leads straightforwardly to the 
following leading-order results: 






D + D l 



D + D* u ) 



D* + D°) 



B(B+ -»■ D* + D*°) 



JD+ 
f D + 


sin 2 9 C B(B+ - 


- DtD°) , 


Id+ 

f 2 


sin 2 6 C B(B+ - 


-> D+D*°) 


9 2 D *+ 

g 2 D *+ 


sm 2 e c B{B+- 


-> z^+d ) 


9 2 D *+ 
9 2 D *+ 
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(Al) 



where #c is the Cabibbo angle, fx and g x » (X = D + or D+) are the decay constants. Since 
our present knowledge of fn+, fo+i 9d*+ an d 9 D *+ is quite poor ||, we take fn+ ~ fo+ ~ 0-8 
and g D *+ ~ 9 n *+ ~ 0.8 f° r simplicity and illustration M. Choosing the central values of 

" s 

B{B+ -► D+Z) ), etc §, we approximately obtain B{B+ -► D+Z) ) « 5.3 x 10~ 4 , 23(5+ -► 
L>+Z>°) « 3.1 x 10" 4 , B{B+ -»■ L>*+L>°) « 3.7 x KT 4 and 8(5+ -> D* + D*°) « 7.1 x 10" 4 . 
From this rough estimation one can see that the above decay modes are definitely detectable 
in the first-round experiments of a i?-meson factory. 

To roughly estimate the CP asymmetry between Z? + — ► D + D° and Z?~ — > D~D°, we 
take the time-like penguin contribution into account [K|. The annihilation and space-like 
penguin effects are expected to be negligible if we insist on the significant formfactor sup- 
pression associated with them []. Then the overall decay amplitudes can be calculated, by 
use of the QCD-improved effective weak Hamiltonian and factorization approximation, in a 



renormalization-scheme independent way [ITJ, [L2|]. Instead of repeating the technical details 
of such a treatment, here we only write out the resultant expressions of Sf and S{ in the 
assumptions made above: 
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(A2) 



^However, one should keep in mind that such an argument may not be on a solid ground and has to be 
examined after some theoretical (experimental) progress is made in deeper understanding of the dynamics 
of nonleptonic B decays. 



where the common hadronic matrix element (D + \(cd)v-A\0) (D°\(bc)v-A\B^) has been sin- 
gled out from 5" and S%. In Eq. (A2), a s and a e are the strong and electroweak coupling 
constants respectively; c, stands for the renormalization-scheme independent Wilson coeffi- 
cient; £ c = 2w|, + /[m c (m& — m c )] arises from the transformation of (V — A)(V + A) currents 
into (V — A)(V — A) ones for Qs-s! and F q (k 2 ) denotes the penguin loop-integral function 
with momentum transfer k at the scale \x = 0(rrib): 



FJk 2 ) = 4 / dx x(l - x) In 
Jo 



m 2 q — k 2 x(l — x) 



mi 



(A3) 



The absorptive part of F q (k 2 ) emerges for k 2 > 4m 2 , leading to the possibility of direct CP 
violation (1^ . 



One can calculate S? and Sf for the decay modes B+ -» D*+D°, D + D*° and D*+D*° 
using the same factorization approximation. If the polarizations of final-state vector mesons 
are summed over, we arrive at the same formulas as Eq. (A2) with £ c = for D* + D°, 
£ c = — 2m 2 D+ /[m c (m b + m c )] for D + D*°, and £ c = for D* + D*°. Of course, such results 
depend upon the assumptions made above and cannot be taken too seriously. 

With the help of Eqs. (A2) and (A3), one is able to evaluate the CP asymmetry A±o 
defined in Eq. (15) and the phase parameter (p given in Eq. (23). For illustration, we 
typically choose m u = 5 MeV, m c = 1.35 GeV, mj, = 5 GeV and m t = 174 GeV. The 
Wolfenstein parameters are taken to be A = 0.22, A = 0.81, p = 0.05 and 77 = 0.36. We 



adopt values of the Wilson coefficients Cj obtained in Ref. [01|. The unknown penguin 
momentum transfer k 2 is treated as a free parameter changing from O.Olm^ to m 2 . A few 
points can be drawn from the explicit numerical calculations: 

(a) The QCD (gluonic) penguin plays the dominant role in the overall penguin amplitude, 
while the electroweak penguin effect is negligibly small. At k 2 = Am 2 w 0.3m| both A±o 
and tan(2<^) undergo a remarkable change in magnitude. 

(b) The CP asymmetries „4 ±0 between B+ -► D^+D^° and B~ -► D^~D^° have 
the same sign and are of the order 3%. The relative change of each asymmetry due to the 
uncertain penguin momentum transfer k 2 is less than 15%. 

(c) With the inputs listed above, the phase parameter tp is estimated to be around —3°. 
Considering the large uncertainties associated with the inputs and the approach itself, we 
believe that a significant deviation of if from zero (e.g., <p ~ —10°) cannot be excluded even 
within the standard model. 

(d) Observation of the above CP-violating signals to three standard deviations needs 
about 10 8 £?„ events, if the composite detection efficiency is at the 10% level. Such mea- 
surements are possible in the first-round experiments at the forthcoming B factories. 



In the case that the decay channels Bd — ► £)(*) £)(*) are significantly suppressed, we ex- 
pect that the direct CP asymmetries in Bd — »■ £)(*)+Z)M~ are comparable in magnitude with 
those in 5+ -»■ D^ + D^° vs £" -»■ D^~D^° (see Eq. (19) for illustration). Nevertheless, 
much more B^B^ events are needed to detect the former on the T(4S f ) resonance due to the 
cost for flavor tagging. It is likely to measure direct CP-violating signals in Bd — ► £)W+Z)(*)~ 
in the second-round experiments of a 5-meson factory. 
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